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ABSTRACT 

- - I Context. The radius of an exoplanet may be affected by various factors, including irradiation received from the host 
(N. star, the mass of the planet and its heavy element content. A significant number of transiting exoplanets have now 

' been discovered for which the mass, radius, semi-major axis, host star metallicity and stellar effective temperature are 

known. 

Aims. We use multivariate regression models to determine the power-law dependence of planetary radius on planetary 
^ ■ equilibrium temperature Teq, planetary mass Mp, stellar metallicity [Fe/H], orbital semi-major axis a, and tidal heating 

■ rate Hudai, for 119 transiting planets in three distinct mass regimes. 

_ Methods. We fit models initially to all 119 planets, resulting in fairly high scatter between fitted and observed radii, 

and subsequently to three subsets of these planets: Saturn-mass planets, Jupiter-mass planets, and high-mass planets. 
Results. We find models for each subset that fit the observed planetary radii well and show the importance of the various 
(N- environmental parameters on each subset. 

Conclusions. We determine that heating leads to larger planet radii, as expected, increasing mass leads to increased 

or decreased radii of low-mass (< 0.5Rj) and high-mass (> 2.0Rj) planets, respectively (with no mass effect on 
Jupiter-mass planets), and increased host-star metallicity leads to smaller planetary radii, indicating a relationship 
between host-star metallicity and planet heavy element content. For Saturn-mass planets, a good fit to the radii may be 
rH ; obtained from log{Rp/Rj) = -0.077 + 0.450 log(Mp/Mj) - 0.314 [Fe/H] + 0.671 log(a/AU) + 0.398 log(re,/K). The 

p radii of Jupiter-mass planets may be fit by \og{Rp/R,j) = -2.217 + 0.856 \og{Tsq/K) + 0.291 log{a/AU). High-mass 

|! planets' radu are best fit by log{Rp/Rj) = -1.067 + 0.380 log{nq/K) - 0.093 log(Mp/Mj) - 0.057 [Fe/H] + 0.019 

O ' log{Htidai/i X lO'^'^). These equations produce a very good fit to the observed radii, with a mean absolute difference 

^ I between fitted and observed radius of 0.11 R,j, compared to the mean reported uncertainty in observed radius of 0.07 R,j. 

C/3 , A clear distinction is seen between the core-dominated Saturn-mass (0.1-0.5 M,/) planets, whose radii are determined 

■ almost exclusively by their mass and heavy element content, and the gaseous envelope-dominated Jupiter-mass (0.5- 
2.0 Mj) planets, whose radii increase strongly with irradiating flux, partially offset by a power-law dependence on 

, I orbital separation. 

^ ' Key words. Stars: planetary systems 

a^: 

ON- 

. 1. Introduction Additional heating effect s have also been pro posed, 

including Ohmic h eating (iBatvein fc StevensonI I201GI : 

. „, ,. 1 1 ■ r 11 1 11 iBatvgin et al.l 120111 : iLaughlin et al.l 1201 iD . a mechanism 

Ihe mass-radius relationship oi gravitationallv bound ob- r T , ^i ,• r ,i t r i.i i 

• . , , • ; 1 1 1 ; 01 planetary heating from the coupling of the plane- 

fNl . lects may be approximate d by a polytropic relatio nship , i:- c u j i, ■ a i • i- i, 4-- 

. i n 71^71-n'i/fs-n'i frr*^ ' i, . \ — HI -, „„„ K rni tary maguetic held and atm ospheric nows; kinetic heating 

T-H; givenby iJoc Md ") ( |Burrows fc Lieb _ed[i99i. The dcuillot fc Showmanl [200l where a small amount of flux 

> polytropic index n ranges from for low mass planets r 77 T 7 7 ■ , r j ■ j- i • j 

' 1 r ■ -11;; • n 71 ; o ;^ r from the host star is transformed into kinetic energy and 

• ^ i made of incompressible matter, i.e. R oc M to 3/2 for ;i ;i i • ^.i, i. r 4.1 i I j 

is/> ^ ' then thermaf energy m the atmosphere of the planet; and 

^ : low-mass stars withelectron degeneracy, i.e. i? oc A/ tidal heating (IBodenheimer et al.l[M . [2051lJackson et all 

H . (IChabrier et al. | 120091). Between these two regimes, at ^qq^ ^ue to circularisation of close-in exoplanet orbits. 

- - ■ appro ximately Jupiter-mass objects l|Zapolskv fc Salpeten 

[T969I) . n ~ 1, i.e. R oc Af°, showing that the radii of Jupiter- The metallicity of the system is also thought to affect 
mass planets do not depend on their masses. The variation planetary radu. Increasing the metallicity could lead to an 
in radii of most known exoplanets must therefore arise from increased heavy element abundance in the atmosphere, pro- 
other inffuencing factors. For example, iGuillotI (|2005() de- ducing higher atmospheric opacity which could again retard 
rive a mass-radius relationship for highly irradiated planets, the lo ss of heat from a pla net and thus slow its contrac- 
incorporating heating effects. tions (jBurrows et al.|[2007f ). However, an enriched atmo- 
sphere would also have a higher mean molecular weight. 
Various potential effects on planetary radii have been offset ting partially or totally any red uction in contrac- 
discussed in the literature. The equilibrium temperature of tion (jBurrows et al.|[2007l : lGuillotl(2008[) . Alternatively, the 
the planet, related to the irradiating flux r eceived from the higher metallicity may result in increased heavy elements 
host star, is likely to affec t the radius ([Guillot et al.l[T996t in the core, leading to a decreased planetary radius rel- 
iGuillot fc Showmanll200^ in that planets that are strongly ative to a planet o f the same mass but lower me tallicity 
heated by their host stars may have inflated radii. ( Guillot et al.ll2006[ l. lFortnev fc NettelmannI (|2010r) set this 
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out in their 'basic questions' about planetary structure and 
composition, including whether heavy elements in a planet 
mix into the H-He envelope, or are found in the core, and 
whether a planet's heavy element mass fraction depends on 
that of its p arent star. 

Recently, iLaughlin et al.l ()2011t ) found a correlation be- 
tween the radius anomalies from comparing observed radii 
with model radii based on calculations of coreless giant 
planets tabulated for var ious mass and insolation values 
(jBodenheimer et aI1l2003[) . They found an additional ra- 
dius dependence on planetary temperature of i? oc T^-^, and 
suggested Ohmic heating may be largely responsible. They 
also noted that even including this additional dependence, 
there remained significant scatter in the observed radii and 
suggested evidence of a signficant correlation with host star 
metallicitics, and that tidal heating may also contribute. 

The system age may also have an effect since the ra- 
dius of an isolated planet should decrease with time due to 
Kelvin-Helmholt z contraction. 

Ad ditionally, [Valencia et al.l (|201Cl[ ) and iJackson et al.l 
(|2010f ) recently investigated CoRoT-7b, concluding that it 
is likely to have initially had a much larger mass than cur- 
rently, possibly including a substantial gaseous component 
which would have given it a larger radius than the 0.15 Rj 
observed today. Close-in planets are susceptible to evapo- 
ration due to intense radiation from the star. A large loss 
of mass (and radius) can occur due to enhanced atmo- 
spheric escape for exoplanets that are very close to their 
stars. Molecules are able to escape even if below the nor- 
mal escape velocity, needing only to reach the Roche Lobe. 
The atmosphere of H-He planets is only loosely bound to 
the planetary core, so evaporation can be quicker than 
the contraction of a planet, which may result in loss of 
the whole envelope for a planet of up to IMj initially 
([Valencia et al.l[201Clt) . A planet which lost a significant frac- 
tion of mass may become stable when the EU V flux from 
the young host s tar dr ops ijValencia et al.|[20To[) . Models by 
iHubbard et al.l (|2007| ) show significant in itial mass loss in 
the first 0.1 Gyr of a svstem. fBaraffe et al.l (|2005il show that 
recent Neptunian objects could have originally been more 
massive irradiated gas giants. 

The reduction in re quired escape energ y may be given 
by (jJackson et al.ll2010t lEAaev et al.ll2007f) 
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where 
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(2) 



The mass loss rate due to stellar irradiation is then given 

by 



dMp 
dt 



(3) 



where F^uv is the extreme UV flux from the star and e 
is the fraction of incoming energy that is removed by es- 
caping gas. Thus planets of smaller semi-major axis may 
be trimmed down by atmospheric blow-off. Atmospheric 
escape of HD 209458, with Kude = 0.65, has been ob- 
served in the Ly a line ([Vidal-Madiar et al]l2003|) . More re- 
cently, iFossati et al.l (|2010[ ) report finding enhanced transit 
depths when observing the very close-in exoplanet WASP- 
12b with the HST/COS instrument using UV transmission 



spectroscopy. All the planets considered here have Kude 
values below 1, while some have very low values and thus 
significant enhancement of mass loss, for example WASP- 
19b and WASP-12b with Kude values of 0.29 and 0.27 re- 
spectively. 

We previously investigated the effect of metallicity and 
equilibrium temperature on the radii of planets in a mass 
range of 0.1 — O.GMj, finding a strong correlation of stellar 
radius with irradiating flu x, and a weaker an ticorrelation 
with host-star metallicity (jEnoch et al.ll20l"ol) . In this pa- 
per, we seek to illuminate the main trends seen in the radii 
of the majority of known transiting exoplanets with a wide 
range of masses and radii using several parameters. Section 
|5]sets out the method of analysis and gives the resulting cal- 
ibrations for the planetary radii, while Section |3| provides a 
discussion of the results. 



2. Analysis 

There are currently (November 2011) almost 200 known 
transiting exoplanetfQ, for which values have been pub- 
lished for both mass and radius. To perform our analysis we 
selected all known transiting planets in the mass range 0.1 - 
12.0 Mj with orbital periods below 10 days for which values 
have been published for host star metallicity and effective 
temperature. This resulted in a sample of 119 planets, given 
in Tables m to H 

We used upda t ed pa rameters for WASP-17b, as given 
in I Anderson et aP (1201 ih . with a planetary radius of 2.0 Rj 
and eccentricity 0.03 foun d from secondary eclipse timing, 
for WASP-7b, as given in ISouthworth et al.l ()20lH) . with a 
planetary radius of 1.33 R,j, and fo r CoRoT-13b, with a 
radius of 1.25 Rj (jSouthworth|[20Tll) . 

We flxed the eccentricities of many of these plan- 
ets to zero: WASP-6b, WASP-12b, WASP-22b, WASP- 
28b and W ASP-32b were rep orted to have eccentrici- 



ties o f 0.05 (iGillon et al.ll2009h . 0. 02 (iHebb et al.ll2009h . 
0.02 (iMaxted et al.l I2010bir a05 (|Westl 120 111) "and 0.02 
(iMaxted et al.ll2010al) respec tively. However, Lucy-Sweeney 
tests (iLucv fc Sweenevlll97lD on the original data give val- 
ues of 0.22, 0.17, 0.26, 1.0 and 0.68 respectively for the prob- 
ability that the measured eccentricities could have arisen 
by chance from underlying circular orbits, given the uncer- 
tainties and uneven sampling of the radial-velocity obser- 
vations, showing that in all these cases a circular orbit is 
likely. 

Secondary eclipse timings often show eccentricities to be 
considerably lower than the upper limits allo wed by radial- 
veloci ty me asurements alone. Fo r WASP-12b. lHusnoo et al.l 
(1201 iD and lCampo et al.l ()2010D find a likely circular orbit, 
based on secondary eclipse timings. Si milarly, WASP- 19b 
was attributed an eccentricity of 0.02 in lHebb et al.l ()2010() 
but iGibson et al.l (|2010[ ) find a circular orbit from the tim- 
ing of the measured secondary eclipse, though conversely, 
a small eccentricity (0.009) in the o r bit of WASP- 18b was 
confirmed to 7a by iNvmever et al.l ()2011l ) from their sec- 
ondary eclipse measurements. WASP-lOb was originally 
repor ted t o have an eccentr i city o f 0.06 ([Christian et al.l 
120091 ) ■ but lMacieiewski et al.l (|2010f ) find that a circular or- 
bit is likely, with the difference in radial velocities due to 
starspot activity. 
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Additionally, iPont et al] (j201ll) point out that a finite 
best-fit eccentricity will always be found where the eccen- 
tricity is allowed to float as a free parameter in fitting the 
data. Spurious eccentricities with apparent significance as 
great as 3a often arise from a combination of measurement 
error and uneven sampling around t he orbit. From a h omo- 
geneous re-analysis of available data. lPont et al.l (j201l|) find 
that the eccentricities of CoRoT-5b and HAT-13b should 
be set to zero. For other planets in our sample announced 
with small eccentricity values, we also fix them to zero if 
the values are within 3a of 0, which was the case for HAT- 
18b, -19b, -20b, -22b, -23b, -25b, -28b, -29b, -32b and - 
33b, CoRoT-12b, -18b and -19b, Kepler-12b, -14b and -17b, 
WASP-23b and -50b and HD 209458b, for which the data 
has p reviously been noted to be consistent wit h a circular 
orbit (|Laughlin et al.ll2005l: [Snellen et al.ll20Q8D . 

This resulted in 17 out of 119 planets having an eccen- 
tricity greater than zero, as shown in Tables [1] to 01 with 
significantly more high-mass planets (greater than 2.0 Mj) 
showing a non-zero eccentricity: 11 of 32 high- mass plan- 
ets, compared to 1 of 16 Saturn-mass (0.1-0.5 Mj) and 4 of 
71 Jupiter-mass (0.5-2.0 Mj) planets. This difference does 
not seem to be explained by the kn own period-eccentricity 
relationship (jHalbwachs et al.ll2005l for example) , that exo- 
planets on short periods have circular orbits due to tidal cir- 
cularisation while exoplanets on wider orbits have a range 
of eccentricity values, since both low-mass and high-mass 
eccentric and non-eccentric planets cover a range of orbital 
periods. The differ ence may therefore lend suppor t to the 
trend discussed by iRibas fc Miralda-Escudel ()2007() for ra- 
dial velocity-discovered planets that high-mass planets have 
an eccentricity distribution closer to binary stars than low- 
mass planets. It should be noted though that it is easier 
to detect genuine eccentricity in high-mass planets than in 
low-mass ones because the velocity measurement errors are 
small relative to the total RV amplitude. 

To determine the factors that have a significant effect on 
planetary radii, we performed multivariate regression anal- 
yses (Singular Value Decomposition (SVD) fits), weighted 
by the uncertainties on the measurements for planetary 
radii, for a variety of explanatory parameters, including the 
planetary equilibrium temperature, Tgg, in Kelvin, semi- 
major axis, a, in AU, planetary mass. A/, in Jupiter masses, 
stellar metallicity, [Fc/H] and tidal heating rate on the 
planet, Hudai, divided by 1 x 10^°W. The equilibrium tem- 
perature of the planet is given by 



eq 




(4) 



where ^g- is the host star's effective temperature, i?* is 
the stellar radius, a is the semi-major axis, A is the plane- 
tary Bond albedo and F is the fraction of the planet surface 
that is re-radiating flux. Albedo values have no t so far been 
determ ined for the majority of exoplanets, but lRowe et al.l 
(I2OOI) determined a very low albedo for HD 209458b, find- 
ing just 0.04 ± 0.05 for the geometric albedo and inferring 
a Bond albedo of less than 0.12 from this. We set A = 
[F ~ 1) here to calculate the equilibrium temperatures. For 
a given A and F then, Teg is directly proportional to the 
irradiating flux from the star. 

The tidal heating rate is given by 

_ 63 (GA/,)3/2A/,i?5a-i5/2e2 



where Q is the tidal dissipation f actor (taken to be 1 0^) 
and k is the Love number (=0.51) (| Jackson et al.l[2008h . 

To determine the best-fitting SVD model, we used the 
Bayesian Information Criterion, BIC, value where 



BIC = x^ + Np\n{ND) 



(6) 



where Np is the number of param eters used in the fit 
and N ]j is the number of datapoints (ISzvdlowski fc Kurekl 
1200^ . Initially, we used a one-parameter model, and used 
the BIC values to determine the parameter with the 
strongest influence on planetary radii. Next, we tried a two- 
parameter model, using the best parameter from the one- 
parameter fits, and adding each of the other parameters in 
turn. If a two-parameter model was found to have a lower 
BIC value than the best one-parameter model, we selected 
this and attempted a three-parameter model using the best 
two-parameter model plus a further parameter, and so on 
until the best model was found. 



2.1. All planets 

To find the best-fitting model for planetary radius for all 
119 planets over the full mass range of 0.1 to 12.0 Afj, we 
initially subtracted the e xpected cold-body m ass-radius re- 
lationship determined bv lSeaeer et al.l (|2007l ) using interior 
models. Figure |T] shows log(M/Afj) against \ o^(R/Rj) for 
the 11 9 planets used in this analysis, with the lSeaeer et al.l 
(120071) mass-radius relationship for pure hydrogen with a 
solid line. Most of the planets follow the pure H relation- 
ship, or lie above it, while a few (notably CoRoT-8b), per- 
haps strongly core-dominated, lie below it. We subtracted 
the logarithmic mass-radius relationship expected for pure 
H from the planets, leaving the \og{M /Mj)-\og{R/ Rj) 
scatter shown in [2] 



a: 



O 




Fig. 1. Mass-Radius plot of 119 transiting exoplanets show- 
ing t he expected mass-ra dius relationship for pure H (solid 
fine) ISeager et all (|2007l ) 
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(5) 
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Fig. 2. Mass against radius residuals after subtracting the 
iSeaeer et alJ ([2007h mass-radius relationship. 



We then performed an SVD analysis on the radius resid- 
uals, and found a best-fit model where 

log(i?p/i?j) = -1.374 + 0.6121og(reg/A') 
+ 0M4\og{a/AU) 

+ 0M3\og{Hudai/l X 10^0 

- 0.023[Fe/H] (7) 

A table of BIG values summarising the model building pro- 
cess is given in Table [5] 

The use of the logarithm of planetary equilibrium tem- 
perature may incorporate several possible heating effects 
showing a power-law relationship of planetary temperature 
with radius. 

A relationship between planetary radius was not only 
found with planetary equilibrium temperature, as expected, 
but also with the semi-major axis. The relationship with 
temperature contains a separate and opposing dependence 
on semi-major axis than that seen in the direct relationship 
with semi-major axis where the planetary radius tends to 
be smaller for closer-in planets. This tendency can actually 
be seen directly by comparing pairs of planets for which pa- 
rameters other than the semi-major axis are very similar. 
For example WASP- 19b and WASP-4b have very similar 
mass (1.15 and 1.12 Mj, respectively), metallicity (~0) and 
equilibrium temperatures (1995 and 1875 K, respectively) 
but WASP-19b has a radius of l.SliJj at 0.016 AU while 
WASP-4b has 1.42i?j at 0.023 AU. Direct comparisons be- 
tween planets are generally difficult due to the different val- 
ues of the various parameters affecting the radii, and thus 
analyses such as set out here can help to disentangle each 
effect. The direct relationship of semi- major axis on radius 
is seen more clearly when analysing sub-samples of planets, 
below, and is discussed further in those sections. 

A plot of observed planetary radii against fitted radii 
using the above equation is given in Figure O showing sig- 
nificant scatter (mean error in fitted radius of 0.16 Rj), 
with 35 of 119 fitted radii values more than 0.2 Rj from 
the observed values (7 with more than 0.4 Rj error). This 
calibration was not therefore felt to provide a very good 
fit to the data. We then partitioned the 119 planets into 16 
Saturn-mass planets, 0.1 — 0.5Mj, dealt with in Scction l^?^ 



71 Jupiter-mass planets, 0.5 — 2.0Mj, discussed in Section 
[23 and 32 high-mass planets, 2.0- 12.0Mj, in Section [2:31 
For each subset of planet s, we perforrned SV D fits as above, 
without subtracting the ([Seaeer et al.l[2007l ) mass-radius re- 
lationships first. We discuss the significant effects on radii 
for each mass range below, and provide a best-fitting equa- 
tion. 
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Fig. 3. Results of the rad ius calibratio n on al l 119 transit- 
ing exoplanets using the [Seager et al.l ([2007[ ) mass-radius 
relationship and a single four-parameter SVD fit. 



2.2. Saturn-mass planets 

Sixteen known transiting exoplanets fall into the mass 
range 0.1-0.5A/j (including WASP-6b). The upper mass 
cutoff of O.SMj was chosen to represent the point where the 
mass-radius relationship changes slope from a steep corre- 
lation in the low mass planets to an almost flat relationship 
in the Jupiter-mass group, thereby splitting the planets into 
distinct groups with different relationships between radius 
and mass. 

The best-fitting model for this mass range was found 
using the procedure described for the global SVD fit above, 
with the model BIG values provided in Table [51 and was 
found to be a four-parameter model, mainly determined by 
planet mass and host star metallicity: 

log{Rp/Rj) = -0.077 + 0.4501og(Mp/M,7) 

- 0.314[Fe/H] +0.6711og(a/AC/) 

+ 0.3981og(re,//v) (8) 

The equilibrium temperature of the planet only minimally 
affects the planet radius here (though it has a strong ef- 
fect on more massive planets, below). We did not include 
tidal heating in the fit, since only 1 of the 16 planets have 
eccentricities demonstrably greater than zero (WASP-17b). 

Eight planets have radii with residuals greater than 
±0.2 Rj relative to the fitted model: HAT-18b (0.20), HAT- 
19b (0.31) and Kepler-12b (0.26) all have their radii un- 
derestimated, such that their observed radii are greater 
than expected from the model fit. The radii of GoRoT- 
5b, GoRoT-8b, WASP-6b, WASP-13b and WASP-21b are 
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overestimated by 0.32, 0.21, 0.20, 0.26 and 0.47 Rj respec- 
tively. 

Planetary mass As mentioned in Section [TJ the radii of 
gaseous planets do not generally depend on their masses but 
there is an effect on radii for lower-mass planets, made up of 
essentially incompressible matter. For the Saturn-mass sub- 
set of planets, the mass has the strongest effect on planet 
radius. A correlation between mass and adjusted radius, af- 
ter removing the effects of the other parameters, produces a 
coefficient of 0.51, shown in Figure|H with a least-absolute- 
deviation model overplotted. The positive coeSicients in the 
correlation and the SVD model fit shows that for low-mass 
planets, radius increases as mass increases, as expected. 

The fit coefficient for mass, at 0.45, is in fact higher 
than expected from the polytropic relationship with n 1 
for incompressible matter, such that R oc M^/^, or log R oc 
0.33 log M . The additional dependence of radius on mass 
indicates that the density of the planet itself is dependent 
on mass, suggesting that the core-envelo pe ratio increases 
towards lower mass, as recently shown in lMiller &: Fortiievl 
(l20Tll) . 
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Fig. 4. The correlation between observed radius and mass 
for Saturn-mass planets, after removing the effects of the 
other parameters used in the SVD fit. 



Metallicity A correlation between the adjusted radius, af- 
ter removing effects due to other parameters, and the host 
star metallicity results in a coefficient of -0.56, shown in 
Figure [5j The strong negative coefficient, as well as the 
negative coefficient value in the SVD fit, indicates that as 
the metallicity of the host star increases, planetary radius 
decreases. This lends support to the theory that higher sys- 
tem metallicii :y produces larger planetary cores, and thus 
smaller radii (jGuillot et al.ll2006h . Given that a strong cor- 
relation is found at all between host star metallicity and 
planetary radius, this implies that a planet's heavy element 
mass fraction is related to that of its host star. 



Planet-star separation A small effect on BIG value was also 
produced by the star-planet separation, however a corre- 
lation between the semi-major axis and the observed ra- 
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-0.6 -0.4 -0.2 -0.0 0.2 0.4 0.6 
[Fe/H] 

Fig. 5. The correlation between observed radius and metal- 
licity for Saturn-mass planets, after removing the effects of 
the other parameters used in the SVD fit. 



dius, after removing the effects on radius of the other terms 
by subtracting these parameters multiplied by their coeffi- 
cients above, produces a low coefficient value of 0.34, shown 
in Figure El clearly not a strong effect in this mass range. 
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Fig. 6. The correlation between observed radius and 
planet-star distance for Saturn-mass planets, after remov- 
ing the effects of the other parameters used in the SVD 
fit. 



Stellar irradiation The equilibrium temperature of the 
planet also has an effect on the radii. Performing a simple 
correlation between the logarithms of observed radius and 
equilibrium temperature, after adjusting for the other terms 
in the radius fit, gives a weak coefficient of 0.28, shown in 
Figure [71 The positive coefficient shows that planets of a 
higher equilibrium temperature generally have larger radii, 
but for core-dominated Saturn-mass planets, this appears 
to be only a minor effect. Th e apparent contrast in results 
found here in comparison to lEnoch et al.l (|2010[ ) is due to 
lo wering the up p er ma ss cutoff from 0.6Mj to 0.5Mj here. 
In lEnoch et all (|2010D . over half the planets in the small 
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sample analysed were in the mass range 0.5-0.6 Mj, which 
is here part of the larger Jupiter-mass planet group where 
a large effect of equilibrium temperature on radius is seen 
(see below). 



1 .4 fe 

; 0.1-0.5 Mj 
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Fig. 7. The correlation between observed radius and equi- 
librium temperature for Saturn-mass planets, after remov- 
ing the effects of the other parameters used in the SVD 
fit. 



Age All 16 planets have published values for isochrone age. 
To see if the age of a system affects the radius of close-in 
planets, for example due to Kelvin-Helmholtz contraction, 
we tried adding a parameter for age to the model, but it 
produced no improvement in the model, with BIC = 474 
compared to 470 without the age parameter. The expected 
Kelvin-Helmholtz contraction is likely to be mostly offset 
by the continued heating from strong stellar irradiation on 
the close-in planets. 

2.3. Jupiter-mass planets 

This is the largest subset of transiting planets, with 71 hav- 
ing a mass between 0.5 and 2.0 Mj. A fit to all 71 planets 
resulted in a two-parameter fit given by 

\og{Rp/Rj) = -2.217 + 0.8561og(Te,/i^) 
+ 0.2911og(a/AL/) 

(9) 

Again, we did not attempt to fit tidal heating, since only 
4 of the 71 planets have non-zero eccentricity. This model 
produced a mean error of 0.106 R,j in fitted radii, with five 
planets having errors in fitted radii of greater than 0.2i?j 
(up to 0.27 Rj), plus one significant outlier, HAT-32b, un- 
derestimated by 0.63 Rj. HAT-32b is an extremely bloated 
planet with a measured radiu s of over 2 i?./, whic h may 
be in a slightly eccentric orbit (jHartman et al.|[2010[ ). Tidal 
heating on the planet could potentially explain the extra 
infiation above that fitted here. 

The addition of a coefficient for mass produced a small 
improvement of 20.1 in BIC value, but increased the mean 
error in fitted radii to 0.113, so was felt to be insignificant, 
particularly since a correlation of mass with radius adjusted 



for the fitted effects of temperature and semi-major axis 
yielded a very small coefficient of 0.13. Model BIC values 
are given in Table [71 

Planetary temperature is the strongest parameter in 
the fit: the high planetary equilibrium temperatures sus- 
tained due to the strong stellar irradiation on close-in ex- 
oplanets produces a shallow temperature gradient in the 
planet, impeding the interior cooling and thus contrac- 
tion. These planets must have migrated early and not con- 
tracte d as a planet orbiting fur ther from the star would 
have (|Guillot fc ShowmanI [2002h since if they had formed 
and cooled at a larger semi-major axis before inward mi- 
gration, their radii would on ly be marginally increased due 
to a puffed-up atmosphere (jlro et al.ll2005t iBurrows et al.l 
l2000f ) and such a strong effect on radius due to irradiation 
would not be seen. 

Correlations of radius with equilibrium temperature and 
semi-major axis, after removing effects on the radius due to 
the other terms, give coefficients of 0.84 and 0.62, respec- 
tively, shown in Figures \E\ and [HI The positive correlation 
and fit coefficient for the semi-major axis indicate that the 
closer the planet is to the star, the smaller the radius. This 
could be due to the ease of atmospheric escape described 
in Section [1] where the atmospheric escape of closer plan- 
ets is enhanced by the Kude factor given in Equation [TJ As 
mentioned above, WASP-19b has a Kade value of 0.29 at 
0.016 AU, and its measured radius is 0.11 Rj lower than 
that of WASP-4b, of similar mass, equilibrium temperature 
and host star metallicity, orbiting slightly farther from its 
host star. 




0.2 f ^ 

2.8 2.9 3.0 3.1 3.2 3.3 3.4 3.5 



Log(T,q (K)) 

Fig. 8. The correlation between observed radius and equi- 
librium temperature for Jupiter-mass planets, after remov- 
ing the effects of semi-major axis. 

In this mass range, planet mass has no apparent effect 
on planet radius. This is as expected from the relationship 
R oc M^^~"^/^'^~"\ as discussed in Section [1] where the 
polytropic index is n ~ 1 for Jupiter-mass planets. A lack 
of mass-radius effect in this mass range also indicates that, 
in contrast to the Saturn-mass planets discussed above, the 
cores of these planets are either unimportant, or are not 
correlated with planet mass. 

Adding a coefficient for age to the model increases the 
BIC value from 218 to 220, for the 57 planets with published 
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in 



9^. -2.7 




-2.0 -1.8 -1.6 -1.4 
Log(a (AU)) 



1.2 



1.0 



log(i?p/i?j) = -1.067 + 0.3801og(Te,//O 

- 0.0931og(Afp/Afj) - 0.057[Fe/H] 

+ 0.0191og(fft,dai/l X 10^0) 



(10) 



with model BIC values provided in Table [51 The mean error 
in fitted radius is 0.10 with four planets having an error 
in fitted radius greater than 0.2 Rj (WASP-18b, CoRoT- 
2b, CoRoT-17b and CoRoT-23b), though none was greater 
than 0.3 Rj. Adding a coefficient for semi-major axis did 
not improve the BIC. 

The correlation of equilibrium temperature and high- 
mass planet radii gives a coefficient of 0.72, shown in Figure 

m 



Fig. 9. The correlation between observed radius and semi- 
major axis for Jupiter-mass planets, after removing the ef- 
fects of equilibrium temperature. 



isochrone ages. To investigate whether the relationship be- 
tween semi-major axis and radius has any correlation with 
the isochrone age of the system, we altered the model pa- 
rameter from log(a) to log(a)/log(a(7e). This increased the 
model BIC value to 306. If there is no dependence here 
on age, a reduction in radius of close-in exoplanets due to 
blow-off coul d occur early in the sy stem history. The T 
Tauri phase (jjovl 119451 : iHerbid Il952f ) in the evolution of 
low- mass stars (< 2Mq) produces strong XUV radiation 
and intense stellar winds that cause the disa ppearance of 
the gas component of the disc (|IIavashilll98lt) via outflow. 
The Ktide enhanced atmospheric escape due to the XUV 
radiation combined with the strong stellar wind could strip 
these close exoplanet s of their outer layers of a tmosphere, 
reducing their radii. iMurrav-Clav et ahl ()2009() calculated 
that a hot Jupiter at 0.05 AU from its host star could lose 
around 6 x lO^^g s~^ due to the intense XUV from a T-Tauri 
star, equating to around 1% of a Jupiter mass during the 
~100 Myr T Tauri phase; more for a planet orbiting closer 
than 0.05 AU. Thus the relationships of radius with semi- 
major axis and equilibrium temperature could both be due 
to proximity of the planet to the host star, but from en- 
tirely different processes, namely atmospheric blow-off and 
retarded contraction, respectively. However, a small reduc- 
tion in the overall mass of a Jupiter-mass planet is not likely 
to produce a noticeable effect on radius in this mass regime, 
so the radius-semi-major axis relationship may be due to 
another factor, perhaps a migration halting mechanism re- 
sulting in planets with larger cores / higher heavy element 
content orbiting closer to their host stars than otherwise 
similar planets with lower heavy element content. 



2.4. High-mass planets 

Thirty-two transiting exoplanets used in this analysis have 
a mass in the range 2.0-12.0 Mj. Equilibrium temperature, 
planet mass, metallicity and tidal heating were found to be 
important in determining their radii. The SVD fit resulted 



O 




2.8 2.9 3.0 3.1 3.2 3.3 3.4 3.5 

Log(T,q (K)) 

Fig. 10. The correlation between observed radius and equi- 
librium temperature for high-mass planets. 



As expected, the mass term again becomes important 
for more massive planets, greater than a few Jupiter masses, 
where electron degeneracy in the dense cores produces a 
contraction of the planet and thus a smaller radius. The 
correlation of planet mass with radius adjusted for equilib- 
rium temperature and semi-major axis effects is -0.47, as 
seen in Figure [TTJ Planet mass here thus has the opposite 
effect on radius as it does on low-mass planets, where the 
correlation and SVD fit coefficient were positive. 

The addition of a term for tidal heating also has a small 
effect on the fit - 12 of the 32 planets have non-zero eccen- 
tricity. A correlation between tidal heating and adjusted 
radius gives a coefficient of 0.28. 

Metallicity also has an effect on planet radius: a radius- 
mctallicity correlation gives a fairly weak coefficient of - 
0.25, shown in Figure 1131 but the inclusion of a term for 
metallicity does improve the SVD fit. 

The semi-major axis is no longer seen to have the same 
effect on planet radii as was seen for lower-mass planets. 
The higher mass of these planets may act to prevent the 
atmospheric escape that may produce the correlation of 
radius with semi-major axis seen in lower-mass planets. 

Adding age to the model slightly decreases the model 
BIC, from 105 to 96, for the 26 planets that have pub- 
lished isochrone ages, so a small age dependence is seen for 
high-mass planets. This SVD fit produces a small negative 
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Fig. 11. The correlation between observed radius and mass Fig. 13. The correlation between observed radius and 
for high-mass planets, after removing the effects of other metallicity for high-mass planets, after removing the effects 
terms. of the other terms. 



-0.95 
-1 .00 



: 2.0-1 2.0 Mj 









-3-2-10 1 2 
Htidoi 

Fig. 12. The correlation between observed radius and tidal 
heating for high-mass planets of non-zero eccentricity, after 
removing the effects of the other terms. 

coefficient for age of -0.03, showing that radius decreases 
as the planet ages, as would generally be expected due to 
Kelvin-Helmholtz contraction. 

3. Discussion 

We have attempted to determine the factors that have an 
effect on planetary radii, along with the magnitude and di- 
rection of each significant effect. A positive relationship be- 
tween planetary equilibrium temperature and radius is seen 
for Jupiter and high-mass planets, showing that the higher 
the equilibrium temperature, the larger the planet radius, 
though not significantly for the Saturn mass planets. The 
radius dependence found here is not as strong as that re- 
ported recently by iLaughlin et al.l (l201lf ) who found a pro- 
portionality of i?p (X Tg^^ above that expected from theoret- 
ical calculations. Here, the Saturn-mass and low-mass plan- 
ets show around i?p oc while the Jupiter-mass planets 
show i?p oc Tg^^, in total. These lower values of propor- 
tionality suggest that kinetic heating (jGuillot fc ShowmanI 



I2OO2I ), expected to lead to oc T^f (|Laughlin et al.l 
I2011D . could be the major source of providing the addi- 
tional heating to inflate planetary radii, since Ohmic heat- 
ing should yield a higher dependence on temperature of 
around i?p oc Tl;^ (|Laughlin et al.ll20lH) . 

A negative relationship between host star metallicity 
and planetary radius is seen for the Saturn-mass and high- 
mass planets, showing that as host star metallicity in- 
creases, planet radius decreases. This implies firstly that 
there is a relationship between host star metallicity and 
the fraction of heavy elements present in a planet, and 
secondly that the higher metallicity likely results in larger 
pla netary cores, prodti cing smaller radii, both as discussed 
by ([Guillot et al.|[2006t ). No relationship with metallicity is 
seen for the gaseous envelope-dominated Jupiter mass plan- 
ets. 

A strong positive relationship between semi-major axis 
and planetary radius is seen for the Jupiter-mass planets, 
such that planets closer to their host stars have smaller 
radii, an effect in opposition to the bloating caused by 
strong stellar irradiation which tends to make planets closer 
to their stars larger. This may be due to atmospheric blow- 
off in the early stages of the system development or a mi- 
gration stopping mechanism that leads to planets of higher 
heavy element contect orbiting closer to the star. The re- 
lationship is much weaker for the core-dominated Saturn- 
mass planets, and no such relationship is seen for the high- 
mass planets, perhaps because their larger gravitational 
pull prevents atmospheric escape. 

The radii of the planets rises with mass for low-mass 
planets, is unaffected by mass for Jupiter-mass planets, and 
falls with mass for the high-mass planets. This changing re- 
lationship is as expected with the move from incompressible 
matter to partially electron degenerate bodies at high mass, 
though the relationship with mass for Saturn-mass planets 
is stronger than expected from the polytropic relationship, 
indicating an additional dependence that may be due to 
the core-envelope ratio increasing towards lower mass. 

A clear contrast between the core-dominated Saturn- 
mass planets of < 0.5Mj and the gaseous envelope- 
dominated Jupiter-mass planets of 0.5-2.0 Mj is therefore 
seen in the fits. Planetary mass and heavy element content 
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almost exclusively determine the radius of a Saturn-mass 
planets, with stellar irradiation having little effect, while 
stellar irradiation and semi-major axis determine the ra- 
dius of a Jupiter-mass planet almost entirely. 

Once we had determined the best SVD model for each 
subset of planets, we performed a Monte Carlo SVD anal- 
ysis of 50,000 runs to determine the uncertainties on the 
fit coefficients. In each run, each input parameter to the 
SVD fit was perturbed randomly on a normal distribution 
about their observed values with standard deviations as 
their published parameter uncertainty values. A summary 
of the radius calibration terms and uncertainties is provided 
in Table H 

Using the fits obtained for the three subsets of transit- 
ing planets to obtain fitted radii for the 119 planets results 
in Figure 1141 with a generally good agreement between fit- 
ted and observed radii. The mean error in fitted radius 
to the observed radius is 0.11 Rj compared to the mean 
uncertainty in reported observed radius of 0.07 Rj. This 
solution using three subsets of planets provides a better 
calibration to radius than the single global analysis, where 
35 planets had fitted radius errors of greater than 0.2 Rj 
(17 greater than 0.3 Rj), and the mean error in fitted radius 
was 0.16 Rj. Here, 18 of the 119 planets have errors greater 
than 0.2 Rj (only 4 greater than 0.3 Rj): 8 Saturn-mass 
planets, 6 Jupiter-mass planets and 4 high-mass planets. 
The two significant outliers are WASP-21b, with an overes- 
timated fitted radius of 1.54 Rj compared to its measured 
radius of only 1.07 Rj, and HAT-32b with an extremely 
large measured radius of 2.04 Rj compared to the low fit- 
ted radius of 1.41 Rj. The simiarly bloated Saturn-mass 
planet WASP-17b, with measured radius 1.99 Rj is fitted 
fairly weU to 1.87 




Fitted radius 

Fig. 14. Results of the radius calibration on 119 transiting 
exoplanets using the 3 equations given for each subset. 

Approximations such as setting the albedo of all plan- 
ets to zero in calculating planetary equilibrium temperature 
may have contributed to the poorer radii fits for some plan- 
ets, and the model scatter may be reduced if all albedos 
were known. Additionally, the planetary metallicity may 
sometimes be enhanced or depleted compared to an aver- 
age planet orbiting a host star of a certain metallicity, which 



may lead to larger or smaller than expected cores and thus 
smaller or larger than expected radii. 
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Table 1. Saturn-mass transiting planets used in the analysis of radii. 



ID 


Mp (Mj) 




Rp (Rj) 


a (AU) 


e 


Teff (K) 


[Fe/H] 




R* (Rq) 


Age (Gy) 


HAT- 18b 


0.20 ±0.01 


1, 


.00 ± 0.05 


0.0559 ± 0.0007 





4803 ± 80 


0.10 ±0.08 


0, 


,75 ± 0.04 


12.4 


HAT- 12b 


0.21 ±0.01 




96+0 03 

'-'•^"-0.02 


0.0384 ± 0.0003 





4650 ± 60 


-0.29 ±0.05 




70+0 02 

"^-0.01 


2.0 


CoRoT-8b 


0.22 ± 0.03 


0, 


.57 ±0.02 


0.063 ± 0.001 





5080 ± 80 


0.3 ±0.1 


0, 


.77 ±0.02 


2.0 


WASP-29b 


0.25 ±0.02 


0, 


.74 ± 0.06 


0.0456 ± 0.0006 





4800 ± 150 


0.11 ±0.14 


0, 


,85 ± 0.05 


5.0 


WASP-39b 


0.28 ±0.03 


1, 


.27 ±0.04 


0.0486 ± 0.0005 





5400 ± 150 


-0.12 ±0.1 


0, 


,90 ± 0.02 


9.0 


HAT- 19b 


0.29 ±0.02 


1, 


.13 ±0.07 


0.0466 ± 0.0008 





4990 ± 130 


0.23 ±0.08 


0, 


,82 ±0.05 


8.8 


WASP-21b 


0.30 ±0.01 


1, 


.07 ±0.05 


0.0520 ± 0.0004 





5800 ± 100 


-0.4 ±0.1 


1, 


,06 ± 0.04 


5.0 


HD149026b 


0.37 ±0.01 


0, 


.81 ±0.03 


0431+0 0007 

U.U40i_g 0006 





6147 ±50 


0.36 ± 0.05 




1 t:4+0.05 


0.2 


Kepler-7b 


0.43 ± 0.04 


1, 


.48 ± 0.05 


622+0 0011 





5933 ± 44 


0.11 ±0.03 


1, 


,84 ± 0.07 


1.0 


Kepler- 12b 


0.43 ± 0.04 


1, 


.70 ± 0.03 


0.0556 ± 0.0007 





5947 ± 100 


0.07 ±0.04 


1, 


,48 ± 0.03 


4.0 


WASP- lib 


0.46 ± 0.03 




1 nc,+o-05 

l.UO_Q Q3 


0.043 ± 0.002 





4980 ± 60 


0.13 ±0.08 


0, 


,81 ±0.03 


4.1 


WASP-13b 


46+0 06 

"•™-0.05 


1, 


.21 ±0.14 


5 27+0-0017 

U.UOZ 1 _o.0019 





5826 ± 100 


0.0 ±0.2 


1, 


,34 ±0.13 


4.9 


CoRoT-5b 


47+0-07 

-0.02 




1 ■30+004 

±.0»_Q Q5 


0.0495 ± 0.0003 





6100 ± 65 


-0.25 ± 0.06 


1, 


,19 ±0.04 


1.4 


WASP-31b 


0.48 ± 0.03 


1, 


.54 ± 0.06 


0.0466 ± 0.0003 





6200 ± 100 


-0.19 ±0.09 


1, 


,24 ± 0.04 


4.0 


WASP-17b 


0.49 ± 0.03 


1, 


.99 ± 0.08 


0.0515 ± 0.0034 


0.03 


6650 ± 80 


-0.19 ±0.09 


1, 


,57 ± 0.06 


2.7 


WASP-6b 




1, 


.22 ±0.05 


0421+0 0008 

U.U4Zi_Q Q0J3 





5450 ± 100 


-0.20 ± 0.09 




r, 07+0.03 


7.0 
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Table 2. Jupiter-mass transiting planets used in the analysis of radii (part 1). 



ID 




A/p {Mj) 




Rp (Rj) 


a (AU) 


e 


Teff (K) 


[Fe/H] 


R* (Rq) 


Age (Gy) 


HAT-lb 


0, 


.52 ±0.03 


1 


.23 ±0.06 


0.0553 ±0.0014 





6047 ± 56 


0.12 ± 0.05 


1.12 ± 0.05 


1.0 


OGLE- 111b 


0, 


.53 ±0.11 


1 


.07 ±0.05 


0.047 ±0.001 





5070 ± 400 


0.12 ±0.28 


0.83 ± 0.03 


- 


HAT- 17b 


0, 


.53 ±0.05 


1 


.01 ±0.05 


0.0882 ±0.0010 


0.35 


5246 ± 100 


0.0 ±0.1 


0.84 ±0.05 


7.8 


WASP-15b 


0, 


.54 ±0.05 


1, 


.43 ± 0.08 


0.0499 ±0.0018 





6300 ± 100 


-0.17 ± 0.11 


1.48 ±0.07 


1.3 


CoRoT-16b 


0, 


.54 ±0.09 


1 


.17 ±0.15 


0.0618 ±0.0015 


0.33 


5650 ± 100 


0.19 ± 0.06 


1.19 ±0.14 


6.73 


WASP-22b 


0, 


.56 ± 0.02 


1 


.12 ± 0.04 


0.0468 ± 0.0004 





6000 ± 100 


-0.05 ± 0.10 


1.13 ± 0.03 


- 


XO-2b 


0, 


.57 ± 0.06 





.97 ± 0.03 


0.0369 ± 0.002 





5340 ± 32 


0.45 ± 0.02 


0.96 ± 0.02 


0.7 


HAT-25b 


0, 


.57 ±0.05 


1 


.19 ±0.05 


0.047 ± 0.001 





5500 ± 100 


0.31 ±0.10 


0.96 ± 0.05 


3.2 


WASP-25b 


0, 


.58 ±0.04 




1 99+0.06 
-'-•^^-0.05 


0.0473 ± 0.0004 





5703 ± 100 


-0.07 ±0.10 


0.92 ± 0.04 


0.1 


WASP-34b 


0, 


.59 ± 0.01 




1 99+0.11 
-'-•^^-0.08 


0.0524 ± 0.0004 


0.04 


5700 ± 100 


-0.02 ± 0.1 


0.93 ± 0.12 


6.7 


HAT-3b 


0, 


.60 ± 0.03 





.89 ± 0.05 


0.0389 ± 0.0007 





5185 ± 46 


0.27 ± 0.04 


0.82 ± 0.04 


0.3 


Kepler-8b 




0.60t°;i^ 


1 


.42 ± 0.06 


rv r\ ^ oo-l-O 0006 

0.0483l:°;^^ 





6213 ± 150 


-0.06 ± 0.03 


1.49 ± 0.06 


1.5 


TrES-lb 


0, 


.61 ± 0.06 


1 


.08 ± 0.03 


0.0393 ± 0.0007 





5250 ± 200 


0.0 ± 0.2 


0.82 ± 0.02 


1.4 


OGLE- 10b 


0, 


.63 ±0.14 


1 


.26 ±0.07 


0.0416 ±0.0001 





5800 ± 100 


0.0 ±0.2 


1.16 ±0.06 


- 


HAT-28b 


0, 


.63 ± 0.04 


1 01 o+0. 11 


0.0434 ± 0.0007 





5680 ± 90 


0.12 ± 0.08 




6.1 


HAT-27/WASP-40b 


0, 


.66 ± 0.03 




1 06+°''^ 

^^•'-'"-0.04 


0.0403 ± 0.0005 





5300 ± 90 


0.29 ± 0.10 


^ /-./-,-4-0 O'l 


4.4 


Kepler- 15b 


0, 


.66 ± 0.09 





.96 ±0.07 


0.0571 ± 0.0009 





5595 ± 120 


0.36 ± 0.07 


0.99 ± 0.07 


3.7 


Kepler-6b 


0, 


.67 ± 0.03 


1 


.32 ±0.03 


n i~. A ,-■ -,-4-0 OOOf) 

0.0457^;!:;;°* 





5647 ± 44 


0.34 ± 0.04 


1 QQ^*^-"^ 
l.Oy_0.03 


1.0 


HAT-4b 


0, 


.68 ± 0.04 


1 


.27 ±0.05 


0.0446 ± 0.0012 





5860 ± 80 


0.24 ± 0.08 


1.59 ± 0.07 


0.6 


HAT-24b 


0, 


.69 ±0.05 


1 


.24 ± 0.05 


0.0465 ± 0.001 


0.07 


6373 ± 100 


-0.16 ± 0.10 


1.32 ±0.05 


2.8 


HD209458b 


0, 


.69 ± 0.02 




1 09+0.02 

^"J^-O.OS 


0.0471 ± 0.0005 





6000 ± 50 


0.00 ± 0.02 


1.15 ± 0.06 


2.0 


HAT-30/WASP-51b 


0, 


.71 ±0.03 


1 


.34 ±0.07 


0.0419 ± 0.0005 





6250 ± 100 


-0.08 ±0.08 


1.33 ±0.03 


1.1 


CoRoT-4b 


0, 


.72 ± 0.08 




1 iq+006 

^■^^-0.05 


0.090 ±0.001 





6190 ± 60 


0.00 ±0.15 


..^4-001 
l-l'J-0.03 


0.3 


WASP-35b 


0, 


.72 ± 0.06 


1 


.32 ± 0.03 


0.0432 ± 0.0003 





6050 ± 100 


-0.15 ±0.09 


1.09 ±0.02 


5.0 


HAT-33b 


0, 


.76 ± 0.12 


1 


.83 ± 0.29 


0.0503 ± 0.0011 





6401 ± 88 


0.05 ± 0.08 


1.78 ± 0.28 


2.4 


HAT-9b 


0, 


.78 ± 0.09 


1 


.40 ± 0.06 


0.053 ± 0.002 





6350 ± 150 


0.12 ± 0.20 


1.32 ± 0.07 


1.4 


HAT-29b 








1 11+014 


0.0667 ± 0.0008 





6087 ± 88 


0.21 ± 0.08 


1-^^-0.08 


2.2 


TrES-4b 


0, 


.84 ± 0.07 


1 


.67 ±0.06 


0.0509 ± 0.0007 





6100 ± 150 


0.00 ± 0.09 


1.74 ± 0.09 


2.0 


HAT- 13b 


0, 


.85 ± 0.04 


1 


.28 ± 0.08 


n n Ai^^4-Ci OOOfi 





5638 ± 90 


0.43 ± 0.08 


1.56 ± 0.08 


0.8 


WASP- 16b 


0, 


.86 ± 0.06 


1 


.01 ±0.07 


n ^ A^-, 4-0 001 

o.o42il:°{;°^« 





5550 ± 130 


0.01 ± 0.10 


0.95 ± 0.05 


2.2 


WASP-23b 


0, 


.88 ±0.10 





.96 ± 0.06 


U.UO/U_Q 0024 





5150 ± 100 


-0.05 ±0.13 


„ ___Ln O"? 
U-' '-0.05 


- 


WASP- lb 


0, 


.89 ± 0.20 


1 


.36 ±0.10 


0.0382 ± 0.0013 





6200 ± 200 


0.26 ±0.03 


1.38 ±0.10 


1.0 


WASP-44b 


0, 


.89 ± 0.06 


1 


.14 ±0.11 


0.0347 ± 0.0004 





5410 ± 150 


0.06 ±0.10 


0.93 ± 0.07 


0.9 


XO-lb 


0, 


.90 ± 0.07 


1 


.18 ±0.04 


0.0488 ± 0.0005 





5750 ± 13 


0.02 ± 0.04 


0.93 ±0.03 


2.0 


WASP-28b 


0, 


.91 ±0.06 


1 


.12 ± 0.06 


0.046 ± 0.0005 





6100 ± 150 


-0.29 ±0.10 


1.05 ± 0.06 


2.0 


WASP-2b 


0, 


.91 ±0.09 


1 


.02 ± 0.01 


0.0307 ±0.0011 





5200 ± 200 


-0.08 ±0.08 


0.83 ±0.08 


- 


WASP-41b 


0, 


.92 ±0.07 


1 


.21 ± 0.07 


0.0400 ± 0.0005 





5450 ± 150 


-0.08 ±0.09 


1.01 ± 0.26 


1.8 


CoRoT-12b 


0, 


.92 ± 0.07 


1 


.44 ±0.13 


0.0402 ± 0.0009 





5675 ± 80 


0.16 ± 0.10 


1.12 ± 0.09 


6.3 


HAT-32b 


0, 


.94 ± 0.17 


2 


.04 ±0.10 







6001 ± 88 


-0.16 ± 0.08 


1.39 ± 0.07 


3.8 


WASP-7b 


0, 


.96 ± 0.13 


1 


.33 ± 0.09 


0.0617 ± 0.0010 





6400 ± 100 


0.0 ± 0.1 


1.43 ± 0.09 


2.4 


WASP-48b 


0, 


.98 ± 0.09 


1 


.67 ±0.08 


0.0344 ± 0.0004 





5990 ± 90 


-0.12 ± 0.12 


1.75 ± 0.09 


7.9 


OGLE-182b 


1, 


.01 ±0.15 




1 10+0. 13 
l-l'J-0.08 


0.051 ±0.001 





5924 ± 64 


0.37 ±0.08 


1 14+0-23 
l-l*-0.06 




WASP-45b 


1, 


.01 ±0.05 




1-1"-0.14 


0.0405 ± 0.0009 





5140 ± 200 


0.36 ±0.12 


n 05+0.09 


1.4 


WASP-26b 


1, 


.02 ± 0.03 


1 


.32 ± 0.08 


0.040 ± 0.003 





5950 ± 100 


-0.02 ±0.09 


1.34 ±0.06 


2.0 


CoRoT-lb 


1, 


.03 ±0.12 


1 


.49 ± 0.08 


0.0254 ± 0.0004 





5950 ± 150 


-0.30 ±0.25 


1.11 ±0.05 




0GLE-211b 


1, 


.03 ± 0.20 




i-'J'J-o.og 


0.051 ±0.001 





6325 ± 91 


0.11 ±0.10 


l-O*-0.07 




WASP-24b 


1, 


.03 ± 0.04 




i-i'-'-o.os 


0.0359 ± 0.0003 





6075 ± 100 


0.07 ±0.10 


1 1 c+0.04 
1-10-0.05 


1.6 
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Table 3. (continued) Jupiter-mass transiting planets used in the analysis of radii. 



ID 




Mp (Mj) 




Rp (Rj) 




a (AU) 


e 


Teff (K) 


[Fe/H] 




R* (Rq) 


Age (Gy) 


HAT-6b 


1, 


.06 ±0.12 


1, 


.33 ± 0.06 


0, 


.0524 ± 0.0009 





6570 ± 80 


-0.13 ±0.08 


1, 


.46 ± 0.06 


0.7 


HAT-5b 


1, 


.06 ±0.11 


1, 


.26 ±0.05 


0, 


.0408 ± 0.0008 





5960 ± 100 


0.24 ±0.15 


1 


.17 ±0.05 


1.8 


XO-5b 


1, 


.08 ± 0.04 


1, 


.09 ± 0.06 


0, 


.0487 ± 0.0006 





5510 ± 44 


0.18 ±0.03 


1 


.06 ±0.05 


0.8 


Qatar- lb 


1, 


.09 ± 0.08 


1, 


.16 ±0.05 


0, 


.0234 ± 0.0003 





4861 ± 125 


0.2 ± 0.1 





.82 ± 0.03 


4.0 


CoRoT-19b 


1, 


.11 ±0.06 


1, 


.45 ±0.05 


0, 


,0518 ± 0.0008 





6090 ± 70 


-0.02 ± 0.1 


1 


.65 ± 0.04 


5.0 


WASP-4b 


1, 


.12 ± 0.09 




^■^■^-0.07 




0.023 ±0.001 





5500 ± 150 


0.0 ±0.2 


1 


.15 ± 0.28 


- 


HD189733b 


1, 


.13 ±0.03 


1, 


.14 ±0.03 


0, 


.0310 ± 0.0006 





4980 ± 200 


-0.03 ±0.04 





.79 ±0.05 




OGLE- 132b 


1, 


.14 ±0.12 


1, 


.18 ±0.07 


0, 


.0306 ± 0.0008 





6210 ± 59 


0.37 ±0.07 


1 


.34 ± 0.08 


_ 


WASP-19b 


1, 


.15 ±0.08 


1, 


.31 ±0.06 




0164+0 0005 

U.U10'±_Q gggg 





5500 ± 100 


0.02 ±0.09 






4.5 


TrES-2b 


1, 


.20 ± 0.05 


1, 


.27 ± 0.04 


0, 


.0356 ± 0.0008 





5850 ± 50 


-0.15 ±0.10 


1 


.00 ± 0.04 


2.7 


OGLE-56b 


1, 


.29 ± 0.12 


1, 


.30 ±0.05 


0, 


,0225 ± 0.0004 





6119 ±62 


0.25 ±0.08 


1 


.32 ± 0.06 


- 


CoRoT-13b 


1, 


.31 ±0.07 


1, 


.25 ±0.08 




0.051 ± 0.003 





5945 ± 90 


0.01 ±0.07 


1 


.01 ±0.03 


1.6 


OGLE- 113b 


1, 


.32 ± 0.19 


1, 


.09 ± 0.03 


0, 


,0229 ± 0.0002 





4804 ± 106 


0.14 ±0.02 





.77 ±0.02 


- 


WASP- 12b 


1, 


.41 ± 0.10 


1, 


.79 ± 0.09 


0, 


,0229 ± 0.0008 





6300 ± 150 


0.3 ±0.1 


1 


.57 ±0.07 


0.8 


WASP-50b 


1, 


.47 ± 0.09 


1, 


.15 ±0.05 


0, 


,0295 ± 0.0009 





5400 ± 100 


-0.12 ±0.08 





.84 ± 0.03 


7.0 


HAT-8b 




1 co+O.lS 






0. 


,0487 ± 0.0026 





6200 ± 80 


0.01 ±0.08 






1.0 


WASP-5b 


1, 


.64 ± 0.08 


1, 


.17 ±0.06 


0, 


.0273 ± 0.0006 





5880 ± 150 


0.0 ±0.2 


1, 


.08 ± 0.04 


1.4 


XO-4b 


1, 


.72 ± 0.20 


1, 


.34 ±0.05 


0, 


.0555 ± 0.0011 





5700 ± 70 


-0.04 ±0.03 


1, 


.55 ± 0.05 


0.6 


WASP-3b 




J^- 'O-0.14 




1 01+0.05 
J^-'JJ^-0.12 




0.0317iroof 





6400 ± 100 


0.0 ±0.2 




qi+0.06 

^■'^^-0.12 


1.4 


TrES-5b 


1, 


.78 ± 0.06 


1, 


.21 ±0.02 


0, 


,0245 ± 0.0001 





5171 ±36 


0.20 ±0.08 





.87 ±0.01 


7.38 


WASP-43b 


1, 


.78 ±0.10 






0, 


,0142 ± 0.0004 





4400 ± 200 


-0.05 ±0.17 





.60 ± 0.04 


0.4 


HAT-7b 


1, 


.80 ± 0.06 






0. 


,0379 ± 0.0004 





6350 ± 80 


0.26 ±0.08 




1-841°;?? 


1.0 


WASP-37b 


1, 


.80 ±0.17 






0, 


.0446 ± 0.0019 





5800 ± 150 


-0.40 ±0.12 


1 


.00 ±0.05 


11 


TrES-3b 


1, 


.92 ± 0.23 


1, 


.30 ± 0.08 


0, 


.0226 ± 0.0013 





5720 ± 150 


-0.19 ±0.08 
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Table 4. High-mass transiting planets used in the analysis of radii. 



ID 


(Mj) 




Rp (Rj) 


a (AU) 


e 


Teff (K) 


[Fe/H] 


R* (Rq) 


Age (Gy) 


HAT-23b 


2.09 ±0.11 


1 


.37 ±0.09 


0.0232 ± 0.0002 





5905 ± 80 


0.15 ± 0.04 


1.20 ±0.04 


4.0 


WASP-46b 


2.10 ±0.07 


1 


.31 ±0.05 


0.0245 ± 0.0003 





5620 ± 160 


-0.37 ±0.13 


0.92 ±0.03 


1.4 


Kepler-5b 


2.11 ±0.06 


1 


.43 ± 0.05 


0.0506 ± 0.0007 





6297 ± 60 


0.04 ± 0.06 


1.79 ±0.05 


0.6 


HAT-22b 


2.15 ±0.06 


1 


.08 ±0.06 


0.0414 ± 0.0005 





5302 ± 80 


0.24 ± 0.08 


1.04 ±0.044 


12.4 


HAT-31b 


^•J- ' 0.08 




1 070-24 


0.0550 ± 0.0150 


0.25 


6065 ± 100 


0.15 ±0.08 




3.17 


WASP-8b 


2.23 ±0.17 




1 1 7+0.18 
^■^'-0.06 


0.0793 ± 0.003 


0.31 


5600 ± 80 


0.17 ±0.07 


0.951°°^ 


1.0 


HAT- 14b 


2.23 ±0.06 


1 


.15 ±0.05 


0.0606 ± 0.0007 


0.10 


6600 ± 90 


0.11 ±0.08 


1.47 ±0.05 


0.4 


WASP-36b 


2.28 ±0.07 


1 


.27 ±0.03 


0.0262 ± 0.0003 





5881 ± 138 


-0.31 ±0.12 


0.94 ±0.02 


3.0 


CoRoT-llb 


2.33 ±0.34 


1 


.43 ± 0.03 


0.044 ± 0.005 





6440 ± 120 


-0.03 ±0.08 


1.37 ±0.03 


2.0 


CoRoT-17b 


2.45 ±0.16 


1 


.02 ± 0.07 


0.0461 ± 0.0008 





5740 ± 80 


0.0 ±0.1 


1.59 ±0.07 


10.7 


Kepler- 17b 


2.45 ±0.01 


1 


.31 ± 0.02 


0.0259 ± 0.0004 





5781 ± 85 


0.26 ±0.10 


1.05 ±0.03 


1.78 


Qatar-2b 


2.49 ± 0.09 


1 


.14 ±0.04 


0.0215 ±0.0004 





4645 ± 50 


0.0 ±0.1 


0.71 ±0.02 


- 


CoRoT-21b 


2.53 ±0.37 


1 


.30 ±0.14 


0.0417 ±0.0011 





6200 ± 100 


0.0 ±0.1 


1.95 ±0.21 


4.1 


WASP-38b 


2.71 ±0.07 


1 


.08 ±0.05 


U.U<00_Q QQQg 


0.03 


6150 ± 80 


-0.12 ±0.07 


1 07+0.05 

-"-"J' -0.04 


- 


CoRoT-23b 


2.80 ±0.25 


1 


.05 ±0.13 


0.0477 ± 0.0038 


0.16 


5900 ± 100 


0.05 ±0.10 


1.61 ±0.18 


7.2 


CoRoT-6b 


2.96 ± 0.34 


1 


.17 ±0.04 


0.0855 ±0.0015 





6090 ± 70 


-0.2 ±0.1 


1.03 ±0.03 


- 


WASP- 10b 


3.06l[5i? 


1 


.08 ± 0.02 


U.UO /i_o. 0013 





4675 ± 100 


0.03 ±0.2 


0.78 ±0.04 


0.2 


CoRoT-2b 


3.31 ±0.16 


1, 


.47 ± 0.03 


0.0281 ± 0.0009 





5625 ± 120 


0.0 ±0.1 


0.90 ±0.02 


- 


CoRoT-lSb 


3.47 ±0.38 


1 


.31 ±0.18 


0.0295 ±0.0016 





5440 ± 100 


-0.1 ±0.1 


1.00 ±0.13 


0.6 


WASP-33b 


3.50 ± 0.6 


1, 


.50 ± 0.02 


0.0256 ± 0.0002 





7400 ± 200 


0.1 ±0.2 


1.44 ±0.03 


- 


WASP-32b 


3.60 ± 0.07 


1, 


.18 ±0.07 


0.0394 ± 0.0003 





6100 ± 100 


-0.13 ±0.10 


1.11 ±0.05 


- 


HAT-21b 


4.06 ±0.16 


1 


.02 ±0.09 


0.0494 ± 0.0007 


0.23 


5588 ± 80 


0.01 ±0.08 


1.11 ±0.08 


10.2 


HAT- 16b 


4.19 ±0.09 


1 


.29 ±0.07 


0.0413 ± 0.0004 


0.04 


6158 ± 80 


0.17 ±0.08 


1.24 ±0.05 


2.0 


CoRoT-20b 


4.24 ±0.23 





.84 ± 0.04 


0.0902 ± 0.0021 


0.56 


5880 ± 90 


0.14 ±0.12 


1.02 ±0.05 


0.1 


OGLE2-L9b 


4.5 ± 1.0 


1 


.61 ±0.04 


0.0308 ± 0.0005 





6933 ± 58 


-0.05 ±0.2 


1.53 ±0.04 


0.3 


HAT-20b 


7.25 ±0.19 





.87 ±0.03 


0.0361 ±0.0005 





4595 ± 80 


0.35 ± 0.08 


0.69 ±0.21 


6.7 


CoRoT-14b 


7.6 ±0.6 


1 


.09 ±0.07 


0.0270 ± 0.0020 





6035 ± 100 


0.05 ±0.15 


1.21 ±0.08 


0.6 


WASP- 14b 


7 70+0.43 
' -'^-0.67 




1 26+°-"** 

J^-^"-0.06 


0.037±SJ;'»^ 


0.09 


6475 ± 100 


0.0 ±0.2 


1.30 ±0.05 


0.25 


Kepler-14b 


8.40 ±0.19 




1 TaO-07 


0.0819g:S 





6395 ± 60 


0.12 ±0.06 




2.2 


HAT-2b 


9.09 ±0.24 




1 16+°"'' 
J^-J^"-0.09 


0.0688 ± 0.0007 


0.52 


6290 ± 60 


0.14 ±0.08 


1 64+°"^ 

^^•"^-0.08 


0.5 


WASP- 18b 


10.43 ± 0.40 


1, 


.17 ±0.08 


0.0205 ± 0.0004 





6400 ± 100 


0.0 ± 0.09 


1.23 ±0.05 


0.53 


X0-3b 


11.79 ±0.59 


1, 


.22 ± 0.07 


0.0454 ± 0.0008 


0.26 


6429 ± 100 


-0.177 ±0.08 


1.38 ±0.08 


0.82 



Table 5. BIC values for SVD fits on all 119 transiting exoplanets. 



# parameters 


log(Te,) 


log(a) 


[Fe/H] 


l0g{H tidal) 


1 


2789 


4669 


4614 


4735 


2: Teg + 




2132 


2796 


2791 


3: Teq + a+ 






2103 


2070 


4: Teg ± a+ [Fe/H] ± 






2042 
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Table 6. BIC values for SVD fits on 16 low-mass transiting exoplanets. 



7^ parameters 


iog(re,) 


log(a) 


log(M) 


[Fe/H] 


1: 


1396 


2355 


900 


1551 


2: M + 


925 


895 




612 


3: M + [Fe/H] + 


566 


531 






4: M + [Fe/H] + a + 


470 









Table 7. BIC values for SVD fits on 71 Jupiter-mass transiting exoplanets. 



# parameters 


iog(re,) 


log(a) 


log(M) 


[Fe/H] 


1 


797 


1287 


1628 


1436 


2: Te, + 




402 


765 


406 


3; Te, + a 4- 






382 


406 



Table 8. BIC values for SVD fits on 32 high-mass transiting exoplanets. 



# parameters log(reg) log(a) log(M) [Fe/H] \og{Htidai) 



1 




239 679 


911 


931 


864 


2: 


Teg + 


243 


201 


222 


241 


3: 


Teq + M + 


205 




192 


193 


4: 


Te, + M + [Fc/H] + 


196 






185 


5: 


Te, + M+ [Fe/H] + Hudai + 


187 









Table 9. Summary of radius calibration coefhcients. 



Mass range (Mj) 


Constant 


log(Te,) 


log(a) 


log(Mp) 


[Fe/HJ 


\0^{Htidal) 


0.1-0.5 


-0.077 ±0.698 


0.398 ±0.201 


0.671 ±0.142 


0.450 ±0.153 


-0.314 ±0.048 




0.5-2.0 


-2.217 ±0.834 


0.856 ± 0.284 


0.291 ± 0.126 








2.0-12.0 


-1.067 ±0.117 


0.380 ± 0.036 




-0.093 ± 0.023 


-0.057 ±0.028 


0.019 ± 0.005 



